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Abstract 

In this paper, we present methods for image compression on the basis of eigenvalue de¬ 
composition of normal matrices. The proposed methods are convenient and self-explanatory, 
requiring fewer and easier computations as compared to some existing methods. Through the 
proposed techniques, the image is transformed to the space of normal matrices. Then, the 
properties of spectral decomposition are dealt with to obtain compressed images. Experimen¬ 
tal results are provided to illustrate the validity of the methods. 
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1 Introduction 

Nowadays, digital images and other multimedia files can become very large in size and, therefore, 
occupy a lot of storage space. In addition, owing to their size, it takes more time to move them 
from place to place and a larger bandwidth to download and upload them on the Internet. So, 
digital images may pose problems if we regard the storage space as well as file sharing. To tackle 
this problem, image compression which deals with reducing the size of an image (or any other 
multimedia) file can be used. Image compression actually refers to the reduction of the amount of 
image data (bits) required for representing a digital image without causing any major degradation 
of the image quality. By eliminating redundant data and efficiently optimizing the contents of 
a file image, provided fhat as much basic meaning as possible is preserved, image compression 
techniques, make image files smaller and more feasible to share and store. 

The study of digital image compression has a long history and has received a great deal of 
attention especially with respect to its many important applications. References for theory and 
practice of this method are lISlfTTUTSlfTSll and ifTOl . to name but a few. 
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Image compression, as well as other various fields of digital image processing, benefits from 
the theory of linear algebra as a helpful tool. In particular. Singular Value Decomposition (SVD) 
is one of the most useful tools for image compression II141I17II . 

The matrix A G can be written in the form of USV*, where U G and V G C"^” are 
unitary matrices, i.e., U*U = Im,V*V = where * denotes complex conjugate transpose and /„ is 
nxn identity matrix. The matrix S is an m x n diagonal matrix in such a way that its nonnegative 
entries are ordered in a non-increasing order (see for example. Theorem 7.3.5 of lUl). With respect 
to the influences of singular values of A in compressing an image, and considering fhe imporfanf 
poinf fhaf fhe singular values of A are fhe posifive square roofs of fhe eigenvalues of matrices A*A 
and AA*, the present study concerns itself with the eigenvalue of the normal matrices A + A* and 
A — A* on the purpose of establishing certain techniques for image compression that are efficient, 
lead to desirable results and need fewer calculations. 

In the next section, we briefly presenf some definitions and concepfs abouf normal mafri- 
ces. Secfion |3] consisfs of fwo subsections in which fhe proposed image compression mefhods 
are explained. In Secfion |4j fhe validify rafes of fhe presenfed image compression schemes are 
invesfigafed and compare fheir efficiencies by experimenfal resulfs. 


2 Normal matrices 

In fhis secfion, we review fhe definition and some properties of normal mafrices. See ||4l|7l and 
fhe references mentioned fhere as fhe suggested sources on a series of conditions on normal ma¬ 
trices. In the next section, we will describe the proposed method on the basis of these presented 
properties. 

A matrix M G is called normal if M*M = MM*. Assuming M as an n-square normal ma¬ 

trix, there exists an orthonormal basis of that consists of eigenvectors of M, and M is unitarily 
diagonalizable. That is, let the scalars Ai,..., counted according to multiplicity, be eigenvalues 
of the normal matrix M and let mi ,..., be its corresponding orthonormal eigenvectors. Then, 
the matrix M can be factored as the following: 

n 

M = UAU* = 'Y^XiUiU*, A = diag(Ai,..., A„), 17 = [mi, ... 

;=1 

where the matrix U satisfies UU* = In- Maintaining the generality, assume that eigenvalues are 
ordered in a non-ascending sequence of magnitude, i.e., |Ai | > IA 2 I... > |A„|. 

It is to be noticed that, if all the elements of the matrix M are real, then M* = M^, where M^ 
refers to the transpose of the matrix M. A square matrix M is called symmetric if M = M^ and 
called skew-symmetric if M = —M^. That symmetric and skew-symmetric matrices are normal 
is easy to see. Also, the whole set of the eigenvalues of a real symmetric matrix are real, but all 
the eigenvalues of a real skew-symmetric matrix are purely imaginary. A general square matrix M 
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satisfies M = B + C, for which the symmetric matrix B = (M+M^) /2 is called the symmetric part 
of M and, similarly, the skew-symmetric matrix C = (M—) /2 is called the skew-symmetric part 
of M. As a consequence, every square matrix may be written as the sum of two normal matrices: 
a symmetric matrix and a skew-symmetric one. We specially use this point in the proposed image 
compression techniques. 


3 Image compression methods 

This section consists of two subsections where methods for image compression are presented using 
normal matrices. To this purpose, the matrix representing the image is transformed into the space 
of normal matrices. Next, the properties of its eigenvalue decomposition are utilized, and some 
less significant image data are deleted. Finally, by returning to the original space, the compressed 
image can be constructed. 

Let A be an « X n matrix to represent the image. Two distinct methods are taken into account. 
First, the symmetric parts of X are dealt with to establish an image compression scheme. This 
procedure can be performed in the same way for the skew-symmetric parts of the matrix X. Next, 
another technique is explained using both symmetric and skew-symmetric parts of the matrix X. 
What is noticeable is that finding the eigenvalues and eigenvectors of a matrix requires fewer cal¬ 
culations than finding its singular values and singular vectors. Moreover, it is possible to calculate 
the eigenvalues and eigenvectors of a normal (especially symmetric or skew-symmetric) matrix by 
explicit formulas and, therefore, may yet again need less computation |l2l[3ll3. 

3.1 Image compression method using the symmetric part of X 

In this subsection, a technique of image compression comes into focus on the basis of the eigen¬ 
value decomposition of the symmetric parts of the matrix X. This can be performed in the same 
way for the skew-symmetric parts of A. It is to be noted that, for the same image, the results 
obtained by these two techniques (i.e. using symmetric or skew-symmetric part) may be different. 

Assume Bx as the symmetric part of the matrix A. The normal matrix Bx can be factored as 
in the following: 

n 

Bx = Ubx^BxUbx = Y,^Bx,i^Bx,i‘^Bx,n ^Bx = diag(AB^ j,..., J, 

1=1 

Now, bearing in mind that the eigenvalues are sequenced in a non-ascending order of mag¬ 
nitude, compress the symmetric part of the image by wiping off the small enough eigenvalues of 
Bx- If k of the larger eigenvalues remains, then there is 

k 

Bx = Y^ hx,,UBxXBx ^ ( 1 ) 

1=1 
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where the total storage for Bx is k{n +1). Here comes the stage of reconstructing the compressed 
image X from its symmetric part. To this purpose, all the elements located above the main di¬ 
agonal of the matrix X are needed, and this calls for n{n — l)/2 storage spaces. Because X can 
be considered as an acceptable approximation of X, let us take all the elements above the main 
diagonal of X as the elements of X. Obviously, the elements located below the main diagonal of 
the matrix X^ should be determined too. In addition, it is inferred from the fact 2Bx = X that 
the elements located below the main diagonal of X can be obtained by subtracting the elements 
below the main diagonal of 2Bx from the elements of X^ . See the following equation where ’V” 
and ” X ” denote the given and unknown entries, respectively. 



X / 


X X 

2Bx = 

X X 

X 

+ 

/ X 

X'^ 


It is clear that the main diagonal elements of X are the same as those on the main diagonal 
of Bx- It is also to be noted that, by this procedure, only the elements located below the main 
diagonal of X are modified, and its other elements are remain untouched. 

Moreover, to reconstruct the compressed image X, the elements located below the main di¬ 
agonal of X may be reserved instead of those above the main diagonal, and then a procedure 
similar to what has been performed newly is to be followed. Indeed, X can be partitioned into 
several segments, and those segments may be reserved provided that X is reconstructed. This may 
be useful specially when some segments of the image have more significance or unchanging (or 
uncompressing) some partitions of the image is desirable during the image compression. 


Finally, it should be pointed out that reconstruction of the compressed image X requires 


n-h 1) 4 


n[n—\) 


storage spaces when the symmetric part of X is dealt with. However, if the 


skew-symmetric part of X is concerned, the diagonal elements of the X cannot be obtained from 
Bx and, therefore, they must be reserved. Reconstructing the compressed image X, thus, requires 


n-h 1) 4 


»(»+!) 


Storage spaces when the skew-symmetric part of X is considered for the image 
compression method. In the next subsection, another image compression technique is provided, 
for which reconstructing the compressed image demands fewer storage spaces. 


3.2 Image compression method using both symmetric and skew-symmetric parts 
ofX 

The previous subsection introduced an image compression scheme dealing with almost half of 
an image. By that technique, more than half of the elements of the image remained unchanged. 
This may be of some usages and advantages. Compressed images obtained by the presented 
method have a high quality. However, compressing of just half of the original image may cause 
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the method to lose its reliability as compared to some other image compression schemes. To tackle 
this problem, a new method is presented here about both symmetric and skew-symmetric parts of 
the matrix X in order to compress the image all over. This new image compression technique is 
found to be of a remarkably high reliability. 

As previously mentioned, every matrix equals the sum of its symmetric and skew-symmetric 
parts. To establish the new method, with the definition of Bx, borne in mind, Cx is used to represent 
the skew-symmetric part of X. The matrix Cx can be written as follows: 

n 

Cx = Uc^Ac^Uc^ = =diag(Acx.,,---,^x,J> 

i=l 

With respect to the method described in the previous subsection, for an integer k <n, set 

k 

Cx = Y^ > ( 2 ) 

(=1 

Through ([T]l and Q, the compressed image will be As in the case of reserving 

the matrix Bx, k{n + 1) storage spaces are required for saving the matrix Cx- As a result, the total 
storage requirement for ^ is 2k(n -|- 1). 


4 Experimental results 

In this section, the validity and the influence of the proposed image compression method are ex¬ 
amined. Let us note that the ideas presented in this paper, can readily be used to establish a block- 
based compression scheme. This scheme concerns dividing of an image into non-overlapping 
blocks and compressing of each block. 

The Peak Signal to Noise Ratio (PSNR) is calculated to measure the quality of the compressed 
image. In the case of gray scale images of size M xN, whose pixels are represented with 8 bits, 
PSNR is computed as follows: 

255 ^ 1 9 

PSNR = lOlogf- ; MSB = - ^ij\\ 

cj 

where A,- y and j refer to the elements of the original and the compressed images respectively. In 
the above relationship, MSB stands for the Mean Square Error between the original image and the 
compressed image pixels. In addition. Compression Ratio (CR) may be calculated as an important 
index to evaluate how much of an image is compressed. CR is the amount of bits in the original 
image divided by the amount of bits in the compressed image; that is. 

Original Image Size 

CR =-. 

Compressed Image Size 
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(a) Lena 


(b) Baboon 


(c) Gold hill 


Fig 1: Original Images 


For the sake of simplicity, Method #1 and Method #2 are branded as image compression techniques 
which use the symmetric and the skew-symmetric parts of the matrix representing the image, in 
the image compression schemes, respectively. In addition, let Method #3 denominate the image 
compression technique described in Subsection 13.21 and let Method #4 stand for the image com¬ 
pression method using SVD. Consequently, the following relationships emerge for CR of these 
methods: 


CRMethod#! — 


k(n+l) + ^)’ 


CRMethod#3 — 


2k{n +1) ’ 


and 


CRMethod#2 


CR-Methodi^4 


(k{n + l) + !^)' 

k{2n + 1) 


In the experiments conducted in this study, three 512 x 512 gray scale images were consid¬ 
ered, including (a) Lena, (b) Baboon and (c) Gold hill presented in Fig[T] The PSNR results are 
shown in Tables [T]|2] and |3l for some integer values of k, for images Lena, Baboon and Gold hill, 
respectively. Also, the CR results are given in Table|4]for a 512 x 512 image. The results obtained 
by these techniques are compared to those achieved by Method #4 in our tables. Furthermore, 
Figures |2] and [3] show the compressed image Lena obtained by the proposed techniques as well as 
image compression method using SVD, Method #4, for k = 50 and k = 100, respectively. 
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Method #1 Method #2 



Method #3 Method #4 


Fig 2: Compressed image Lena for k = 50. 



Method #3 Method #4 


Fig 3: Compressed image Lena for k = 100. 
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Table 1: PS NR results for Lena 


k 

Method#! 

Method#2 

Method#3 

Method #4 

10 

21.7012 

22.3391 

22.0050 

22.4065 

20 

24.8722 

25.0305 

24.9462 

25.2937 

30 

26.8832 

27.0315 

26.9531 

27.2243 

40 

28.3856 

28.5746 

28.4751 

28.7927 

50 

29.7358 

29.9001 

29.8129 

30.1761 

75 

32.5608 

32.8031 

32.6763 

33.1093 

100 

35.0504 

35.1668 

35.1047 

35.6641 

150 

39.2993 

39.2955 

39.2938 

39.8988 


Table 2 

: PS NR results for Baboon 

k 

Method#! 

Method#2 

Method#3 

Method #4 

10 

19.7071 

19.7560 

19.7270 

19.8493 

20 

20.6441 

20.7135 

20.6736 

20.7510 

30 

21.3880 

21.4336 

21.4053 

21.5203 

40 

22.0773 

22.1009 

22.0842 

22.2453 

50 

22.7319 

22.7355 

22.7286 

22.9306 

75 

24.2968 

24.3207 

24.3042 

24.5787 

100 

25.8111 

25.7420 

25.7723 

26.1761 

150 

28.7678 

28.6798 

28.7194 

29.3173 


5 Conclusion 

In this paper, the eigenvalue decomposition of symmetric and skew-symmetric matrices, as two 
important kinds of normal matrices, and their properties were applied to obtain image compression 
schemes. The proposed method is straightforward and uncomplicated ones requiring clear and 
fewer computations as compared to some exiting methods. This image compression method is 
the first one ever introduced and concerns the symmetric part of an image, which can also be 
applicable in the case of the skew-symmetric part of the image. The method is capable of keeping 
half (or some selected segments) of an image unchanged, and this feature may be of some usages. 
However, as observed in this study, just half of the original image could be changed, which causes 
the compressed image to lose its reliability. This is what makes the technique different from 
some other image compression techniques. Hence, the second image compression scheme using 
both symmetric and skew-symmetric parts of the original image was proposed. The experimental 










Table 3: PS NR results for Gold hill 


k 

Method#! 

Method#2 

Method#3 

Method #4 

10 

23.9494 

23.9081 

23.9223 

24.1270 

20 

26.3107 

26.2783 

26.2884 

26.6046 

30 

27.8179 

27.8131 

27.8092 

28.1740 

40 

29.0044 

28.9622 

28.9768 

29.3828 

50 

30.0126 

29.9975 

29.9978 

30.4103 

75 

32.1111 

32.1537 

32.1254 

32.6597 

100 

33.9084 

33.8903 

33.8927 

34.5614 

150 

37.0963 

37.0375 

37.0626 

37.8911 


Table 4: CR results for an image of size 512x512 


k 

Method#! 

Method#2 

Method#3 

Method #4 

10 

1.9283 

1.9211 

25.5501 

25.5750 

20 

1.8582 

1.8515 

12.7750 

12.7875 

30 

1.7930 

1.7867 

8.5167 

8.5250 

40 

1.7322 

1.7264 

6.3875 

6.3938 

50 

1.6754 

1.6699 

5.1100 

5.1150 

75 

1.5485 

1.5438 

3.4067 

3.4100 

100 

1.4394 

1.4354 

2.5550 

2.5575 

150 

1.2617 

1.2586 

1.7033 

1.7050 


results show the high reliability of this method. Finally, it is to be pointed out that the proposed 
method can be used to devise block-based compression techniques. 
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